Cross-correlations mediated by Majorana bound states 
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We consider the correlated parallel transport through two quantum dots which are tunnel-coupled to the ends 
of a semiconductor nanowire where the Majorana bound states (MBSs) may emerge under proper conditions. 
In terms of the cross-correlation of currents, we reveal unusual behaviors originated from the nonlocal MBSs, 
including such as the distinct symmetry and antisymmetry of the spectral density in response to the dot-level 
modulations, and the vanished cross correlation occurred when any of the dot-levels is in resonance with the 
Majorana zero mode. 
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It has attracted considerable attention in the past years 
to search for Majorana fermion [[[], ^ in solid states, such 
as in the 5/2 fractional quantum Hall system^] and the 
p-wave superconductor and superfluid[Q ^]. In particu- 
lar, it was predicted that a semiconductor nanowire, with 
strong spin-orbit coupling and subject to external mag- 
netic field, may support zero-energy Majorana bound state 
(MBS) excitation, when the nanowire is in proximity to an 
s-wave superconductor[^-^[|. In this case the hybrid sys- 
tem is driven into a topological superconducting phase. 
The major advantage of this proposal is its involving only 
the most conventional materials, being thus more accessi- 
ble in experiments. Indeed, in recent experiments [ffi-p"2"|] , 
evident signatures of the MBSs were spotted in the hybrid 
nanowire system. 

The realization and manipulation of Majorana fermion 
in solid-states paves a way for the desirable fault-tolerant 
topological quantum computation [[l3|— [l~6|] . In the aspect 
of pure physics, the Majorana study in solid states is of 
interest owing to its many unusual effects, such as the 
non-abelian statistics [ 17 1, the sharp jump in conductance 
peak [jig], the peculiar noise behaviors Jl9|-p^] , and the 
4ir periodic Majorana- Josephson currents |5j |3| p4ft , 
etc. Remarkably, most of these effects are associated 
with the nonlocal nature of the Majorana fermion [|25|-[27|l . 



That is, taking the nanowire realization as an example, the 
emerged pair of MBSs constitute the real and imaginary 
parts of an ordinary fermion. As a consequence, if an elec- 
tron with energy smaller than the energy gap between the 
Majorana zero mode and other exited states is injected into 
the nanowire, the electron would be "split" into two Ma- 
jorana bound states which are essentially correlated while, 
however, spatially separated. 

In this work, we attempt to exploit the Majorana's non- 
locality in a rather direct way. As schematically shown in 
Fig. 1, we consider a setup of parallel transports through 
two quantum dots (QDs), while the two dots are coupled, 
respectively, to two Majorana modes which may appear at 
the ends of the semiconductor nanowire under proper con- 
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FIG. 1 : (Color online) Schematic setup for a cross-correlated 
parallel transport through two quantum dots which are weakly 
tunnel-coupled by a semiconductor nanowire. The nanowire is 
in contact with an s-wave superconductor. Under appropriate 
conditions, such as with large Zeeman splitting and strong spin- 
orbit interaction, a pair of Majorana bound states is anticipated 
to emerge at the ends of the nanowire. Here we explicitly display 
the closed circuit that satisfies the law of current conservation. 



ditions. We will focus on the cross correlation of currents, 
mediated by the nonlocal MBSs, to reveal some unusual 
behaviors. 

Model and Method. — The transport setup is schemati- 
cally shown in Fig. 1, where the semiconductor nanowire 
is in proximity to a conventional s-wave superconduc- 
tor and subject to a Zeeman field. When the Zeeman 
splitting energy V z , the proximity-induced order param- 
eter A, and the chemical potential [i satisfy the condition 
V z > \/ A 2 + /i 2 , the nanowire is predicted to be driven 
into a topological superconducting phase and a pair of 
zero-energy MBSs will emerge at the ends of the nanowire 
[^-|^] . As a result, the two quantum dots which are tunnel- 
coupled to the ends of the nanowire are now effectively 
coupled to the MBSs. The whole system is described by 



H = H 



leads 



Ht + H s , 
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Hamiltonian of the source (S) and drain (D) leads for 
transports through the left (j = 1) and right (j = 2) dots. 



H 



t = £j=i, 2 E a =s,D Efc t ja d}c kja + H.c. describes 
the tunneling between the leads and the dots. The last 
term is the summed system Hamiltonian of the quantum 
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dots and the tunnel-coupled MBSs [[l|]-|20[ ||, 



H sys = e i d ) d i + 9 



£M7i72 



.7 = 1,2 



+ Ai(4 -rfi)7i+«'A 2 (4 + rf2)72- (2) 

In the above Hamiltonians, dj(cij) and ct ^.(cjafc) are, re- 
spectively, the electron creation (annihilation) operators 
of the quantum dots and the attached transport leads, with 
energies Cj and £kj a - 7i and 72 are the Majorana opera- 
tors, while ej\/ ~ e~ e ^ is the overlap coupling amplitude 
between the MBSs, with £ the wire length and £ the super- 
conducting coherence length. 

For the sake of calculation convenience, we move from 
the Majorana to regular fermion (/ and /t) representation, 
via the transformation of -fi = f'+f and 72 = — /). 
The system Hamiltonian is accordingly rewritten as 

H sys = e M Pf + ^ {< ,'/,'/, 

J=l,2 



A J -[(dt+7 7 d j )/ + H.c]}. 



(3) 



Here we introduced a parameter, 77, to distinguish the Ma- 
jorana (r/ = 1) from a regular bound state (RBS, r} = 0). 
This will be useful for latter comparison. Associated with 
the above system Hamiltonian, we introduce the state ba- 
sis as \n1Tifn2), where ni,2,/ = or 1 stand for the oc- 
cupation number of the QDs and the MBSs. One may no- 
tice that the tunneling terms in H sys only conserve charge 
modulo 2e. This reflects the fact that a pair of electrons 
can be extracted out from the superconductor and can 
be also absorbed by the condensate. As a consequence, 
the Hilbert space for the central system can be split into 
two subspaces, {|100), |010), |001), | 111)} with odd par- 
ity and {|011), 1 101) , |110), |000)} with even parity. 

In this work we employ the master equation approach to 
study the cross-correlation between the transports through 
the quantum dots. For a quantum transport setup, the leads 
can be in general regarded as a generalized environment, 
while the central device is the system of interest which is 
described by a reduced density matrix p(t). Moreover, 
under finite bias voltage and in weak transmission limit, 
the Born-Markov approximation applies reasonably, for 
the setup of_Fig. 1 which results in the following master 
equation | 



P 



-ic P + J2 



TfV[dj]p), (4) 
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where Co = [H sys , p] and V[A)p = ApA^ - ±{A+ A, p}. 
In Eq. (Q) we introduced a couple of tunneling rates, 
rj , which are associated with the tunnel-coupling am- 
plitudes and the density-of-states of the leads. 

Based on the solution of the above master equation, the 
currents through each dot can be simply calculated as 



I 



lf = eVfn h 



(5) 



where the occupation of the corresponding quantum dot is 
calculated by rij = Tr[djdj p(t)]. In particular, using the 



stationary density matrix (p), the steady-state currents can 
be obtained. 

Moreover, the master equation approach allows us 
to study the various correlations [|28|], very conveniently. 
For the purpose of the present work, we introduce the 
cross correlation between the currents through the differ- 
ent dots: Gyfah) = Wi^Ijih)], where £[•••] 
stands for the ensemble average in the quantum trajec- 
tory approach[|28||. In this work, we just focus on the 
cross correlation of the drain currents, which are denoted 
here simply by I jit) for brevity. Physically, the cur- 
rent is defined as Ij{t) = e dNj(t)/dt and the aver- 
age is determined by E[dNj(t)] = dtTt[d]dj p{t)]. In 
steady state, the cross correlation is a function of the time 
difference owing to the time-translation invariance, i.e., 
Gij(t2, ti) = Gij(t 2 —ti) = Gij(t). Therefore we have 



Gij(t) = e 2 Tr{rfd4e"(rf^)]dj} 

= e 2 rfrfTr[^Kt]. 



(6) 



Here the Liouvillian propagator e ct describes the state 
evolution governed by Eq. (Q). From the post-selected 
state after a jump, i.e., Pj(0) = djpdj, we introduced 

Pj(t) = e ct pj(0) which obviously satisfies Eq. (Q). 

The spectral density of cross correlation is given by the 
Fourier transform of Gij (t), which yields 

Sijiu) = e 2 rf If Tr[dift(w)4 + dfp^-^d)]. (7) 

In deriving this result, we employed the symmetry prop- 
erty C?2i(i) = Gis{—t), and introduced pj(±uj) = 
Jo Pj(t)e ±lw *<it. Straightforwardly, via the technique of 
Laplace transformation, pj (w) can be solved from Eq. (Q). 
Moreover, in this work we would like to use a sym- 
metrized spectral density, S(u>) = [S^ifw) + S 12(1*})]/ 2, 
to characterize the cross correlation, which reads 

S(u) = e 2 rfrfRe{Tr[d 1 p 2 (o;)4 + d 2 pi(uj)d 2 }}. (8) 



Also, we would like to introduce the Fano factor to char- 
acterize the zero-frequency correlation, which is defined 
as F = S{Q)/2eyfth, where J 1(2) - [Jf (2) + I° (2) )/2 
is the average steady state current through the left (right) 
quantum dot. 

In Refs. [ |26[ |Z/] ], the isolated dot-MBSs-dot system was 
studied. Remarkably, it was found that, even in the limit 
£m — 0, a novel process of teleportationlike electron 
transfer between the distant dots is possible, because of 
the inherent nonlocal nature of the Majorana fermion. If 
€ m 7^ 0, it also leads to a discontinuity of the tunnel- 
ing conductance through one of the dots, under the MBS- 
mediated influence of the other dot. In the following, we 
will show: (i) In the limit cm = 0, there is no cross- 
correlation between the transport currents through the two 
individual dots. This is an interesting result, by noting the 
presence of electron transfer between the dots, mediated 
by the MBSs. (ii) For cm 7^ 0, distinct Majorana signa- 
tures can appear in the cross correlation of currents. 
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FIG. 2: (Color online) Cross-correlation of the steady state cur- 
rents, plotted by the scaled difference of the source and drain 
currents, j3j = (If - If) /If, with j = 1(2) denoting the left 
(right) dot. The cross-correlation is manifested by modulating 
rf/rf and ei of the left dot, via their effects on 02 of the right 
dot. Parameters: Ai = A2 = A, em = £2 = A, rf = rf = A, 
andT? = 0.01A. 
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FIG. 3: (Color online) Transient dynamics of the right dot in 
response to an electron jump from the left dot to its drain lead, 
based on the steady state. The occupation of the jth dot is ob- 
tained by rij(t) = Tr[djdjp(t)], with p(t) evolved from the 
specific "initial" state p(0) = d\pd\/ r Tx(d\pd\). Assuming 
Ai = A2 = A, we give the result in (a) for em = and in 
(b) for 6m = 0.5A. Other parameters: ei = €2 — A, and 

rl% = 0.1 a. 



Stationary and Transient Correlation. — Let us con- 
sider first the steady-state currents through the two dots 
and their correlation mediated by the MBSs. We no- 
tice that, owing to coupling to the MBS, in the transport 
through a quantum dot the source current is in general 
not equal to the drain current [pl|. Their difference, in 
essence, is compensated by a current flowing through the 
nanowire into or from the superconductor condensate. For 
our present double-dot system, owing to the similar rea- 
son, we introduce fa = (If - If) I If (j = 1,2) to 
qualify this unique feature. In Fig. 2, we use this factor 
to illustrate the cross correlation. Specifically, we alter 
rf /rf , and examine its effect on fa and in particular on 
fa. The results in Fig. 2 correspond to cm 7^ 0, where 
the existence of cross correlation is obviously indicated, 
while its peculiar feature is to be further exploited in the 
following. We would like to mention that, if 6m = 0, 
altering rf/rf would have no effect on fa. This result 
completely differs from what happens for two quantum 
dots coupled through a regular bound state (RBS). 

To further illustrate the MBS-mediated cross correla- 



tion, we now turn to the transient dynamics. To be spe- 
cific, we consider how a disturbance on the left dot would 
affect the transport through the right one. Based on the 
steady state, the disturbance can be introduced by a tun- 
neling event of an electron from the left dot to its drain 
lead. This consideration defines a particular "initial" state, 
d 1 pd\/ r T*(dipd\). Then, we consider the transient cur- 
rents associated with this new state. 

First, in Fig. 3(a) we show the result in the limit of 
6m = 0. We find that the right current is unaffected by 
the disturbance on the left dot, which is in agreement with 
the above tunnel-rate modulation. This result straightfor- 
wardly indicates a zero cross-correlation between the cur- 
rents through the two dots. This behavior is indeed not 
obvious at all, by two observations: (i) If we consider the 
two dots being coupled by a regular bound state (RBS), 
even for non-aligned levels of t\ = 62 = A and €m = 0, 
we find that the cross-correlation between I\ and I2 is 
nonzero. ( ii) It was found in Refs. [[HI]], where the iso- 
lated dot-MBSs-dot system was analyzed, that the state 
1 100} can evolve to |001) in short timescale even in the 
limit cm = 0. This seemingly indicates that the right cur- 
rent I2 should be affected by the left-dot disturbance in 
this limit. 

In Fig. 3(b) we show the result of em 7^ 0, assuming 
£a/ = 0.5A and t\ = 62 — A. In the transformed basis, 
this setup does not look too much different from €m = 0. 
However, we find in Fig. 3(b) that the current through the 
right dot is considerably affected by the disturbance on the 
left dot. Moreover, as to be shown in the following, its cor- 
relation behavior differs significantly from that mediated 
by a regular bound state with this same energy (cm). 

Cross Correlation Spectrum. — Below we employ the 
technique of spectral analysis, Eqs. (6)-(8), to analyze the 
transient dynamics as illustrated in Fig. 3(b). In Fig. 4 we 
contour-plot the spectral density of the cross correlation 
in terms of frequency dependence and dot-level modula- 
tions. In Fig. 4(a) and (b) we show the results for the 
MBSs (cm ^ 0), while the results in (c) and (d) corre- 
spond to a RBS with the same energy (cm), as a com- 
parison. We notice that all the results in Fig. 4(a)-(d) are 
an even function of frequency. This is originated from 
the property of Gy (t) = Gji(-t) and the symmetrization 
of the cross-correlation function which leads to Eq. (|8|). 
More interestingly, for the cross correlation mediated by 
MBSs, the spectral density would have different symme- 
try in response to different manner of the dot-level modu- 
lations. That is, corresponding to a single-dot-level mod- 
ulation, the spectral density shown in Fig. 4(a) is antisym- 
metric with respect to the Majorana zero-mode energy, 
while in Fig. 4(b) the spectrum is symmetric under the si- 
multaneous double-dot modulation of the energy levels. 
We find that this is an important signature from the Ma- 
jorana 's nonlocal correlation. In Fig. 4(c) and (d), the 
cross-correlation spectrum, mediated by a regular bound 
state, are always symmetric with respect to the both types 
of dot-level modulations. 

To show more clearly the symmetry behavior, we plot 
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FIG. 4: (Color online) Spectral density of the cross correlation 
versus two types of dot-level modulations: in (a) and (c) modula- 
tion of the left dot level (ei) only; and in (b) and (d) simultaneous 
modulation of the both dot levels (ei = £2 = e). For the single 
dot ei -modulation, we assume £2 = em = A, where A is the 
common coupling strength between the QD and the MBS (i.e., 
Ai = A2 = A). In (a) and (b) we show the results for the MBSs 
with em = A, while the results in (c) and (d) correspond to a 
regular bound state (RBS) with energy em, as a comparison. In 
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FIG. 5: (Color online) Fano factor of the cross correlation, which 
corresponds to the zero-frequency spectrum in Fig. 4. All the 
parameters and explanations are the same as presented there. 



in Fig. 5 the zero-frequency Fano factor. We find that the 



results for MBSs, either symmetric or antisymmetric, are 
centered at the Majorana zero mode and have a vanished 
cross Fano factor there. This feature cannot be interpreted 
by the usual resonance transfer picture in terms of match- 
ing the energy levels. Actually, in the transformed repre- 
sentation, the two MBSs are "combined" to form a reg- 
ular fermion with energy cm- This is likely to give that 
the cross correlation, including the frequency-dependent 
spectrum in Fig. 4(a) and (b), should be centered at the 
"effective" fermion level cm- Also, the vanished cross 
correlation when one of the dot levels (or both levels) is 
in resonance with the Majorana zero mode is far beyond 
our simple intuition. As a sharp comparison, in Fig. 5 we 
also display the results from RBS-mediated correlation. 
Very clearly, we see that the cross Fano factor is symmet- 
ric with respect to the RBS energy for the both types of 
dot-level modulations. Moreover, at the symmetry center 
the cross correlation is nonzero. The remarkable differ- 
ences discussed above, together with the entire different 
"lineshape" in response to the dot-level modulations, pro- 
vide us very useful information to identify the emerged 
Majorana fermion in the nanowire. 

Summary. — We investigated the correlation mediated 
by a pair of Majoranan bound states, by considering par- 
allel transport through two quantum dots. The central part 
of the setup, say, the dot-MBSs-dot system, is experimen- 
tally accessible using the state-of-the-art technology, par- 
ticularly based on a scheme of semiconductor nanowire 
realization. In terms of the cross-correlation of currents, 
we revealed some unusual behaviors, including such as 
the symmetry and antisymmetry of the spectral density in 
response to the dot-level modulations, and the vanished 
cross correlation occurred when any of the dot-levels is in 
resonance with the Majorana zero mode. These important 
signatures, originated from the nonlocality of Majorana 
fermion, are of interest and beyond the simple understand- 
ing based on resonance transfer interpretation. 
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